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Abstract 



In an earlier work made by the first author with J. Turi (Degenerate Dirichlet Problems Related to the 
Invariant Measure of Elasto-Plastic Oscillators, AMO, 2008), the solution of a stochastic variational in- 
equality modeling an elasto-perfectly-plastic oscillator has been studied. The existence and uniqueness of 
an invariant measure have been proven. Nonlocal problems have been introduced in this context. In this 
work, we present a new characterization of the invariant measure. The key finding is the connection between 
■ nonlocal PDEs and local PDEs which can be interpreted with short cycles of the Markov process solution 
, of the stochastic variational inequality. 

Resume 



Une approche analytique de la theorie ergodique des inequations variationnelles stochas- 
tiques. Dans un travail precedent du premier auteur en collaboration avec Janos Turi (Degenerate Dirichlet 
■ Problems Related to the Invariant Measure of Elasto-Plastic Oscillators, AMO, 2008), la solution d'une in- 
\ equation variationnelle stochastique modelisant un oscillateur elastique-parfaitement-plastique a ete etudiee. 
L'existence et l'unicite d'une mesure invariante ont ete prouvees. Des problemes nonlocaux ont ete intro- 
duits dans ce contexte. La conclusion importante est la connexion entre des EDPs nonlocales et des EDPs 
locales qui peuvent etre interpreters comme les cycles courts du processus de Markov solution de l'inequation 
variationnelle stochastique. 



Version francaise abregee 

La dynamique de l'oscillateur elastique-parfaitement-plastique s'exprime a l'aide d'une equation a me- 
moire (voir (l)-(2)). A. Bensoussan et J. Turi ont montre que la relation entre la vitesse et la composante 
elastique de l'oscillateur est un processus de Markov ergodique qui satisfait une inequation variationnelle 
stochastique (voir (3)). La solution admet une mesure invariante caracterisee par dualite a l'aide d'une 
equation aux derivees partielles avec des conditions de bord non- locales (voir (4)). Dans ce travail, une 
nouvelle preuve de la theorie ergodique est presentee ainsi qu'une nouvelle caracterisation de l'unique dis- 
tribution invariante. Dans ce contexte, nous deduisons des nouvelles formules reliant des equations aux 
derivees partielles avec des conditions de bord non-locales a des problemes locaux (voir (10)). 
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1. Introduction 



In the engineering literature, the dynamics of the elastic-perfectly-plastic (EPP) oscillator has been 
formulated as a process x(t) which stands for the displacement of the oscillator, evolving with hysteresis. 
The evolution is defined by the problem 

x + c x + F(x(s),0 < s < t) = to (1) 

with initial conditions of displacement and velocity ir(0) = x, x(0) — y. Here cq > is the viscous damping 
coefficient, k > the stiffness, w is a Wiener process; F(x(s),0 < s < t) is a nonlinear functional which 
depends on the entire trajectory {x(s),0 < s < i) up to time t. The plastic deformation denoted by A(t) 
at time t can be recovered from the pair (x(t), F(x(s), < s < t)) by the following relationship: 

r kY if x(t)=Y + A(t), 

F(x{s),0 < s < t) = I k(x(t)-A(t)) if x(t) e] -Y + A(t),Y + A(t)[, (2) 
[ -kY if x{t) = -Y + A(t). 

where Y is an elasto-plastic bound. Such elasto-plastic oscillator is simple and representative of the elasto- 
plastic behavior of a class of structure dominated by their first mode of vibration, they are employed to 
estimate prediction of failure of mechanical structures. Karnopp & Scharton [4] proposed a separation be- 
tween elastic states and plastic states and introduced a fictitious variable z(t) :— x(t) — A(t). 

Recently, the right mathematical framework of stochastic variational inequalities (SVI) modeling an 
EPP oscillator with noise has been introduced by one of the authors in [2]. Although SVI have been already 
studied in [1] to represent reflection-diffusion processes in convex sets, no connection with random vibration 
had been made so far. The inequality governs the relationship between the velocity y(t) and the variable 

Z (ty. 

dy(t) = -(c y(t) + kz(t))dt + dw(t), {dz(t)-y(t)dt)(<f>-z(t))>0, V|0| < Y, \z(t)\ < Y. (3) 
Let us introduce some notations. 

Notation 1. D := R x (-Y,+Y), D+ := (0, oo) x {Y}, D~ := (-oo,0) x {-Y}, and the differential 
operators A( := -±( yy + (c y + kz)( y -y( z , B + ( := -\C vv + {c y + kY)( y , B_( := -±( yy + (c y-kY)( y . 
where £ is a regular function on D. 

In [2], it has been shown that the probability distribution of (y(t),z(t j) converges to an asymptotic proba- 
bility measure on D U D + U D~ namely v. Moreover, v is the unique invariant distribution of (y(t),z(t)). 
In addition, from [3] we know also that there exists a unique solution u\ to the following partial differential 
equation (PDE) 

Xu x + Au x = f in D, \u x + B + u x = f in D + , \u x + B-u x = f in D~ (4) 

with the nonlocal boundary conditions given by the fact that u x (y, Y) and u x (y, — Y) are continuous, where 
A > and / is a bounded measurable function. The function u x satisfies ||ua||oo < , u x is continuous 
and for all (77, () e D, we have lim^o ^u x (i], () = v(f). We use the notation u x (y 7 z; f). 

Now, we introduce the short cycles to provide a new proof of the ergodic theory for (3). In this context, 
we derive new formulas linking PDEs with nonlocal boundary conditions to local problems. 

1.1. Short cycles 

Let A > 0, consider v x {y, z) the solution of 

\v x +Av x = f in D, Xv x +B+v x = f in D+, \v x + B-v x = f in D~ (5) 
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with the local boundary conditions v\(0 + ,Y) = and v\(0~,— Y) = 0. Also, if / is symmetric (resp. 
antisymmetric) then v\ is symmetric (resp. antisymmetric). We use the notation v\(y,z;f). As A — > 0, 
v\ — > v with 

Av = f in D, B+v = f in D+ , B_v = f in £>" (P v ) 

with the local boundary conditions v(0 + ,Y) = and v(0~ , —Y) = 0. We use the notation v(y,z;f). We 
call v(y, z; f) a short cycle. We detail the solution of (P v ) in the next section. We introduce next Tr + (y, z) 
and 7T~ (y, z) such that 

An + =0 in D, tt+ = 1 in D + , tt + = in D~ (6) 

and 

Att~ = in L>, tt~ = in D + , ir~ = 1 in £>". (7) 

We have 7r + + 7r~ = 1, so the existence and uniqueness of a bounded solution to (6) and (7) are clear. A 
new formulation of the invariant distribution is given by the following theorem. 

Theorem 1.1 (New formulation of the invariant distribution v). Let f be a bounded measurable func- 
tion on D, we have the following analytical characterization of the invariant distribution: 

v(fT,Y;f)+v(0+,-Y;f) 



2w(o+,y ; i) 

_-Y;f) A „ \ , n 

2v x (0-,Y;l) 



Denote u x (f) := ^o-,Y,f)^ -Y if ) ^ ^ A ^ 



u x (y,z:f)-^p-^u(y,z-J), v x {f) v{f) (8) 

where u satisfies 

Au = f-v{f) in D, B+u = f - v{f) in D+ , B_u = f - v{f) in D~ (9) 
with the nonlocal boundary conditions given by the fact that 

u(y,Y) and u(y,—Y) are continuous. 
Then, we obtain also the representation formula 

u(y, z; f) = v(y, z; f) v(f)v(y, z; 1) + ^^l^-'y)^ > Y > V w ( 0+ ' ™ ^ 

2. Analysis of the short cycles 

We describe the solution of (P v ). We can write v(y, z; f) = v e (y, z; f) + v + (y, z; f) + v~(y, z; f) with 
v e ,v + ,v~ satisfying 

Av e = f(y,z) in D, v e = in D + , v e = in D~ , (11) 
Av + = in D, v + (y,Y) =<p + (y;f) in L>+, v + = in D~ , (12) 

and 

Av~=0 in D, v~=0 in D + , v~(y, —Y) = ip~(y; f) in D~ . (13) 
where ip + (y; f) and <fi~(y; f) are defined by 
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^+ y + (c y + kY)<p+ = f(y,Y), y > 0, <^+(0+;/)=0 (14) 

and 

t-2 i oi„vcw pi+y 



-l<Pw + (<%V-kY)<p-=f(y,-Y), y<0, ip~ (0~; f) = 0. (15) 



We check easily the formula <p+(y; f) = 2 d£ exp(-(c £ 2 + 2fcF£)) I? V f(C Y) cxp(-2c £(C - 0)d(, if 
y > and also ^_ (y; f) = 2 / °° ^cxp(-( Co £ 2 - 2kY^)) f(C, -Y) cxp(-2 Co £(C - 0R> if V < 0- 
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2.1. Solution to Problem (11) 

The proof will be based on solving a sequence of Interior Exterior Dirichlet problems and a fixed point 
argument. Thus, we need to state the two following lemmas as preliminary results. It is sufficient to consider 
/ = 1, with no loss of generality. 

2.1.1. Interior Dirichlet problem 

We begin with the interior problem, let Di := (—y~i,y~\) x (—Y,Y), D± := [0, y~\) x {Y}, D± := 
(— yi, 0] x {— Y}. Let us consider the space of continuous functions on [-Y, Y] which are on Y and the 
space Cf of continuous functions on [— Y, Y] which are on —Y. Let ip + 6 and ip~ e Cf . We consider 
the problem 

-^Cyy + (coy + kz)(y - yCz = I in D x , C(y,Y)=0 in D+ , ((y,-Y)=0 in Df (16) 

with ((yi,z) = f + (z) and ((— y~\, z) = <p~(z), if —Y < z < Y. 

Lemma 2.1. There exists a unique bounded solution to the equation (16). 

PROOF. It is sufficient to prove an a priori bound. For that we can assume <p + ,<f~ = 0. Consider A > 
and the function 6(y, z) — exp(Aco(j/ 2 + kz 2 )) then — \Q yy + (c y + kz)6 y — y9 z = 9 (— Ac + 2Acgj/ 2 (l — A)). 
Set next H := —(6 + () then 

-\Hyy + (c y + kz)H y -yH z = -l + 6 (\c - 2Ac 2 l y 2 (l - A)) . (17) 

If we pick A > max(l, ^-) the right hand side of (17) is positive. Therefore the minimum of H can 
occur only on the boundary y = y~\ and z = Y with y > or z = —Y with y < 0. It follows that 
H(y,z) > -cxp(\c Q (yf + Y 2 )) and thus also < C < exp(Ac (y 2 + Y 2 )). 

2.1.2. Exterior Dirichlet problems 

Now, we proceed by considering two exterior Dirichlet problems. Let < y < yi, we define D y<y := 
{y > y, -Y < z < Y}, D± <y := {y > y, z = Y} and D y< _ y := {y < -y, -Y < z < Y}, 
Dy < _y := {y < —y, z = —Y} and consider 

~\vm + ( c oy + kz)r]+ - yrj+ = 1 in D v<y , ri+(y,Y) = in D^ <y (18) 

with the condition r] + (y, z) — ((y, z) if —Y < z <Y, and 

-\%y + (°oy + kz)r]y - yt]~ = 1 in D y< _y, r]-(y,-Y) = in D y< _ y (19) 

with the condition rj~(—y,z) = ((—y,z), if —Y < z < Y. We use the same notation f](y,z) for the two 
problems (18), (19) for the convenience of the reader. We have 

Lemma 2.2. For any y > there exists a unique bounded solution of (18), (19). 

PROOF. It is sufficient to prove the bound, we claim that ||C||oo < < ||CI|oo + ^y^' for y > y and 

IICIU < V(y,z) < IICI|oc + l: ± £ , for y < -y. Consider for instance p(z) = ||C||oo + 2: ^ for y > y, -Y < z < Y 
then -\p vv + py(c y + kz)-yp z = | > 1, p(y,z) = ||CI|i~ + > C(V,z), p(V,z) = IICIU- > 0. So clearly 
V(y> z ) < p( z )- So in all cases we can assert that ||r?||oo < HCIloo + 
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Ve(y,z) 



2.1.3. Solution to Problem (11) 

Proposition 2.1. There exists a unique bounded solution to Problem (11). 

Proof. Uniqueness comes from maximum principle. Setting $ = (ip + (z), ip~ (z)) and using the notation 
= <fi + (z), z) = <p~(z), we can next define T$(y l7 z) = 77(3/1,2) and T$(-yi,z) = r](-yi,z). 

We thus have defined a map F from C^Cf into itself. If T has a fixed point, then it is clear that the 
function 

C(y,z), -yi < y < 3/1, 
v(y,z), y>y, y<~y 

is a solution of (11) since Q = r\ for y < y < 3/1,2 G (— Y,Y) and for — y\ < y < —y,z G (— Y,Y) 
and the required regularity is available at boundary points y, j/i, — y, — j/i. The result will follow from the 
property : T is a contraction mapping. This property will be an easy consequence of the following result. 
Consider the exterior problem 

~ip V y + ip y (coy + kz) - yi> z = in D y<y , ij){y,Y) = Q in D± <y , (20) 

where i/)(y, z) = 1 if — Y < z < Y, then swp_ Y<z< Y ^(Hii z) < p <1. 

Indeed if sup_ y<z<y tp(yi, z) — 1, then the maximum is attained on the line y — y\, and this is impossible 
because it cannot be at z = Y, nor at z = —Y, nor at the interior, by maximum principle considerations. 

2.2. Solution to Problems (12) and (13) 

We now consider the function ip + and ip~ solution of (14) and (15). Note that if y < 0, we have 
(p~(y; 1) = p + (—y; 1). So it is sufficient to consider (14) and we easily see that 



<p + (y;l)=[ cxp(-( Co e 2 + 2fcrO) 1 CX P< 2co ^> dg, if y>0 
Jo ^c £ 

and we have <p + (y; 1) < ^ \og{ °° v ^ Y ), if y > 0. We next want to solve the problem (12). We proceed 
as follow. We extend ip + for y < 0, by a function which is C 2 on R and with compact support on y < 0. It 
is convenient to call ip(y) the C 2 function on M, with compact support for y < and <p(y) = <p + (y; 1) for 
y > 0. We set w + (y, z) = v + (y, z) — ip(y) then we obtain the problem 

Aw + =g in L>, w + (y,Y)=0, in L>+, w + (y, -Y~) = -p(y), in D~ (21) 

with g(y, z) = - (-^Vyy + ( c oy + kz)ip y ). 

But, g{y,z) = l{ 3/ >o} (-1 + k(Y - z)ip y (y)) + l {y<0 } (-(-^f yy + (c Q y + kz)ip y )) and thus, taking into 
account the definition of tp when y < 0, we can assert that g(y, z) is a bounded function. Again, from the 
definition of f(y) when y < 0, we obtain that on the boundary, w + is bounded. It follows from what was 
done for Problem (11) that (21) has a unique solution. So we can state the following proposition. 

Proposition 2.2. There exists a unique solution to (12) of the form v + (y,z) = ip + (y)l{ y>0 y + v + (y,z) 
where v + (y, z) is bounded. Similarly, there exists a unique solution to (13) of the form v~(y, z) = ip~ (y)l^ y<Q y + 
v~(y,z) where v~(y,z) is bounded. 

PROOF. We just define <p{y) extension of (fi + (y) for y < as explained before and consider w + (y, z) solution 
of (21). We know that w + (y,z) is bounded and we have v + (y,z) — ip(y) + w + (y,z) = ip + (y)l{ y>a y + 
ip(y)l{ y<0 } + w + (y, z) which is of the form (12) with v + (y, z) = ip(y)l[ y<0 y + w + (y, z). 

2.3. The complete Problem (P v ) 

Finally, we consider the complete Problem (P v ), we can state 

Theorem 2.1. There exists a unique solution of (P v ) of the form v(y, z; f) = f + (y; f)l{ y> o}+<f~ {y; f)l{ y< o} + 
w(y, z) where w(y, z) is a bounded function which can be written as w — v e + w + + w~ . 

Proof. We just collect the results of Propositions 2.1 and 2.2. 
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3. Ergodic Theorem 

PROOF (Proof of Theorem 1.1). We first prove the result when / is symmetric. In that case, we can 
write 

u x (y, z; f) = v x (y, z; /) + (j v x {y, z; 1)) (22) 

Indeed, we know that u x (y, z; /) and v x (y, z; /) are symmetric. Setting u\(y, z; f) = u\(y, z; /) - v x (y, z; /), 
we obtain 

\u\ + Au x = in D, \u x + B + u x = in D+, \u x + B_u x = in D~ (23) 

with the boundary conditions u A (0+, Y; f)—u x (0~, Y; f) = v x (0~ , Y; f) and u x (0+ , —Y; f)—u x (0~, —Y; f) = 
— v x (0 + , —Y\ /). This last condition is automatically satisfied, thanks to the previous one and the symme- 
try. The function j — v x (y,z; 1) satisfies the three partial differential equations on D, D + and D~ . So, 
u x = C ( j — v x (y, z; 1)) and writing the first boundary condition, we have ii x (0 + ,Y; f) — u x (0~,Y; /) = 

-C{v x (0+,Y;1) - v x (0-,Y;l)) = Cv x (0~ ,Y;1). Hence, C = ^$^rj and formula (22) has been ob- 
tained. Now, we have v x (f) ->■ v{f) = ^q-'p.f) , as A -)• 0. If we define u^(y, z; /) = u\(y, z; f) - ^jp- = 
v x (y,z; })-v x (f)v x (y,z; 1). The function u* x (y,z; f) -> v(y, z; f)-u(f)v(y, z; 1) = v(y, z; f-u(f)), A ->■ 0. 
Also from its definition the function u\(y 7 Y;f) and u^(y, — Y"; /) are continuous. From the choice of v(f) 
the function v(y,Y;f — v{f)) is continuous. Now, since / — v(f) is symmetric v(0 + ,—Y; f — v{f)) — 
v(Q~ , — Y; / — v(f)) = v(0 + , Y; f — v(f)) — v(0~ , Y; f — v(f)) = 0. So the result is proven when / is sym- 
metric. We now consider the situation when / is antisymmetric. We know that u x (y, z; /) is antisymmetric. 
Similarly v x (y,z;f) is antisymmetric. Consider 71T and defined by 

\ir++An+=0 in D, Att+ + B+ir+ = in D+, tt+ = in D~ (24) 

with the boundary condition 7r^(0 + , Y) = 1 and 

Ati\ + AtiT = in D, tt7=0 in £>+, Att~ + B + ir~ = in L>~ (25) 

with the boundary condition 7r7(0~, —Y) = 1. We have 7iT(y, z) = 7T^(— y, — z), we then state the formula 
"A(y, z; /) = v\(y, z; /) ~ ^TI^(^~^^5TT^y^ • So we see that u ^(y, z; /) converges as A -> 0,without 
substracting a number . The function u x (y,z;f ) converges pointwise to u(y,z;f) = v(y,z;f) — 
(n ( y ,z) ^_^y0}^_ ■ y *i) ^ g Q w j len j j s antisymmetric, the results (8)- (9) hold with v x {f) = and u(f) = 0. 
For the general case, we can write / = / sym + / asym with f sym (y,z) = f(y> z )+f(~y- z ) ; f asym (y,z) = 

/(„,«)-/(-„,-,) _ We ^ K/sym) = «(0-,_^) and thus K/sym) = " (0 "' 2 y ;ff + + ;_ ( ^- y;/) Since K/asym) = 0, 

we deduce ^(/) = v{f sym ) = '2v(o+ V -y-i) Y '^ ■ ^ e obtain also the representation formula u(y,z;f) = 
v(y, z; /) - v{f)v{y, z; 1) + n+ t-\o-~Y)' Z) ( V (°~> Y '> f) - w ( 0+ ' - Y '> /)) and the result is obtained. 
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